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ABSTRACT 

One of the main challenges facing standard hierarchical structure formation models is that the pre- 
dicted abundance of galactic subhalos with circular velocities v c ~ 10—30 km s" 1 is an order of magnitude 
higher than the number of satellites actually observed within the Local Group. Using a simple model for 
the formation and evolution of dark halos, based on the extended Press-Schechter formalism and tested 
against N-body results, we show that the theoretical predictions can be reconciled with observations if 
gas accretion in low-mass halos is suppressed after the epoch of reionization. In this picture, the observed 
dwarf satellites correspond to the small fraction of halos that accreted substantial amounts of gas before 
reionization. The photoionization mechanism naturally explains why the discrepancy between predicted 
halos and observed satellites sets in at v c ~ 30 km s , and for reasonable choices of the reionization 
redshift (z IC ~ 5 — 12) the model can reproduce both the amplitude and shape of the observed velocity 
function of galactic satellites. If this explanation is correct, then typical bright galaxy halos contain 
many low-mass dark matter subhalos. These might be detectable through their gravitational lensing 
effects, through their influence on stellar disks, or as dwarf satellites with very high mass-to-light ratios. 
This model also predicts a diffuse stellar component produced by large numbers of tidally disrupted 
dwarfs, perhaps sufficient to account for most of the Milky Way's stellar halo. 

Subject headings: cosmology: theory - galaxies: formation 



1. INTRODUCTION 

Cold Dark Matter (CDM) models with scale-invariant 
initial fluctuation spectra define an elegant and well- 
motivated class of theories with marked success at explain- 
ing most properties of the local and high-redshift universe. 
One of the few perceived problems of such models is the 
apparent overprediction of the number of satellites with 
circular velocities v c ~ 10 — 30 km s _1 within the viri- 
alized dark halos of the Milky Way and M31 (Klypin et 
al. 1999a, hereafter KKVP99; Moore et al. 1999; see also 
Kauffmann et al. 1993 and Gonzales et al. 1998). The 
most recent of these investigations are based on dissipa- 
tionless simulations, so the problem may be rephrased as 
a mismatch between the expected number of dark mat- 
ter subhalos orbiting within the Local Group and the ob- 
served number of satellite galaxies. In order to overcome 
this difficulty, several authors have suggested modifica- 
tions to the standard CDM scenario. These include (1) 
reducing the small-scale power by either appealing to a 
specialized model of inflation with broken scale invariance 
(e.g., Kamionkowski & Liddle 1999) or substituting Warm 
Dark Matter for CDM (e.g., Hogan 1999), and (2) allow- 
ing for strong self-interaction among dark matter parti- 
cles, thereby enhancing satellite destruction within galac- 
tic halos (e.g., Spergel and Steinhardt 1999; but see the 
counter- argument by Miralda-Escude 2000). In this paper, 
we explore a more conservative solution. 

Many authors have pointed out that accretion of gas 
onto low-mass halos and subsequent star formation are 
inefficient in the presence of a strong photoionizing back- 
ground (Ikeuchi 1986; Rees 1986, Babul & Rees 1992; Ef- 
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stathiou 1992; Shapiro, Giroux, & Babul 1994; Thoul & 
Weinberg 1996; Quinn, Katz, & Efstathiou 1996). Moti- 
vated by these results, we investigate whether the abun- 
dance of satellite galaxies can be explained in hierarchical 
models if low-mass galaxy formation is suppressed after 
reionization, and we propose that the observable satel- 
lites correspond to those halos that accreted a substantial 
amount of gas before reionization. This solution is similar 
in some ways to the idea that supernova feedback ejects 
gas from low-mass halos (Dekel & Silk 1986), but it of- 
fers a natural explanation for why some satellite galaxies 
at each circular velocity survive, while most are too dim 
to be seen. In addition, the reionization solution seems 
physically inevitable, while the feedback mechanism may 
be inadequate in all but the smallest halos (Mac Low & 
Ferrara 1998). 

Our approach to calculating the satellite abundance uses 
an extension of the Press-Schecter (1974) formalism to 
predict the mass accretion history of galactic halos and 
a simple model for orbital evolution of the accreted sub- 
structure. For our analysis, we adopt a flat CDM model 
with a non-zero vacuum energy and the following param- 
eters: il m = 0.3, Ha — 0.7, h = 0.7, as — 1.0, where as 
is the rms fluctuation on the scale of Sh" 1 Mpc, h is the 
Hubble constant in units of 100 km s _1 Mpc , and fl m 
and J7a are the density contributions of matter and the 
vacuum respectively, in units of the critical density. 
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Fig. 1. — Average cumulative velocity function of dark 
matter halos accreted by halos of mass 1.1 x lO 12 ^ 1 M 
(z = 0). The average is over 300 merger histories. The thin 
dashed line represents velocity function of all accreted halos; 
the thick dashed line corresponds to halos that survive effects 
of dynamical friction; the thin solid line corresponds to halos 
that survive both the dynamical friction and tidal disruption. 
The thick solid line is a fit to results of cosmological simu- 
lations presented in KKVP99. The figure shows that a sub- 
stantial fraction of accreted substructures can be destroyed by 
z — 0; most of the destruction occurs relatively early (zZ.0.5). 



2. METHOD 
2.1. Modeling the Subhalo Distribution 

Because N-body simulations that resolve low-mass sub- 
halos are computationally expensive, we develop an ap- 
proximate analytic model for the accretion history and 
orbital evolution of satellite halos within a typical Milky 
Way-size dark halo. Our general procedure is to use the ex- 
tended Press-Schechter method (Bond et al. 1991; Lacey 
& Cole 1993, hereafter LC93) to construct the subhalo 
accretion history for each galactic halo and then to deter- 
mine which accreted subhalos are dragged to the center 
due to dynamical friction and which are tidally destroyed. 
The subhalos that survive at z = are used to construct 
the final velocity function of satellite halos. 

We assume that the density profile of each halo is de- 
scribed by the NFW form (Navarro, Frenk, & White 1997): 
Pnfw( x ) k a; _1 (l + x)~ 2 , where x = r/r s , and r s is a 
characteristic inner radius. Given a halo of mass M v ; r at 
rcdshift z, the model of Bullock et al. (1999) supplies the 
typical r s value and specifies the profile completely. The 
circular velocity curve, v 2 (r) = GM(r)/r, peaks at the ra- 
dius r max ~ 2.16r s . Throughout this paper, we use v c to 
refer to the circular velocity at r max . 

We begin by constructing mass growth and halo accre- 
tion histories for an ensemble of galaxy-sized dark mat- 
ter halos. This is done using the merger tree method of 
Somerville & Kolatt (1999, hereafter SK99). We modified 



the method slightly to require that at each time step the 
number of progenitors in the mass range of interest be close 
to the expected average. This modification significantly 
improves the agreement between the generated progeni- 
tor mass function and the analytic prediction. We start 
with halos of mass M vir = 1.1 x 10 12 h^ 1 M e , at z = 0, 
corresponding to v c = 220 km s _1 , and trace satellite ac- 
cretion histories back to z = 10 using time steps chosen 
as described in SK98. We track only accreted halos that 
are more massive than M m = 7 x 10 6 h~ 1 MQ, which corre- 
sponds to v c ~ 10 km s _1 at z = 10 and v c ~ 4 km s _1 at 
z = 0, and treat the accretion of halos smaller than M m 
as diffuse mass growth. 

For each step back in time, the merger tree provides a 
list of progenitors. The most massive progenitor is identi- 
fied with the galactic halo and the rest of the progenitors 
(with M > M m ) are recorded as accreted substructure. 
We are left with a record of the mass growth for the galac- 
tic halo, M v i r (z), as well as the mass of each accreted sub- 
halo, M a , and the redshift of its accretion, z a . For the 
results presented below, we use 300 ensembles of forma- 
tion histories for galactic host halos. Our results do not 
change if the number of ensembles is increased. 

Each subhalo is assigned an initial orbital circularity, 
e, defined as the ratio of the angular momentum of the 
subhalo to that of a circular orbit with the same energy 
e = J/J c . We choose e by drawing a random number uni- 
formly distributed from 0.1 to 1.0, an approximation to the 
distribution found in cosmological simulations (Ghigna et 
al. 1998). Our results are not sensitive to the choice of 
the above range. To determine whether the accreted halo's 
orbit will decay, we use Chandrasekhar's formula to cal- 
culate the decay time, t df , of the orbit's circular radius 
i? c , as outlined in Klypin et al. (1999b). The circular ra- 
dius is defined as the radius of a circular orbit with the 
same energy as the actual orbit. Each subhalo is assumed 
to start at R c = 0.5R v i T (z a ), where R v i r (z a ) is the virial 
radius of the host halo at the time of accretion. This 
choice approximately represents typical binding energies 
of subhalos in simulations. The amplitude of the resulting 
velocity function is somewhat sensitive to this choice. If 
we adopt R c = 0.75i? v ir(2 a ) (less bound orbits), the ampli- 
tude increases by ~ 50%. A change of this magnitude will 
not affect our conclusions, discussed in §4. Once t df is 
known, the fitting formula of Colpi et al. (1999) provides 
the appropriate decay time for the given circularity 

tdf = T% F e QA . (1) 

If tdf is smaller than the time left between z a and z = 0, 
tdf < to — t a , then the subhalo will not survive until the 
present, and it is therefore removed from the list of galactic 
subhalos. 

If tdf is too long for the orbit to have decayed com- 
pletely (tdf > to — t a ), we check whether the subhalo 
would have been tidally disrupted. We assume that the 
halo is disrupted if the tidal radius becomes smaller than 
r max . Such a situation means either a drastic reduction 
in the measured v c (quickly pushing v c below our range of 
interest) or that the halo (and any central baryonic compo- 
nent) should become unstable and rapidly dissolve (Moore, 
Katz & Lake 1996). 

The tidal radius, r t , is determined at the pericenter of 
the orbit at z = 0, where the tides are the strongest. For 
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an orbit that has decayed to a final circular radius R®, the 
pericenter at this time, R p , is given by the smaller of the 
roots of the bound orbit equation (e.g., van den Bosch et 
al. 1999) 



R 



2MR°) - <i>(r)} 
e 2 v 2 {R° c ) 



(2) 



Here, $(r) = —A.6v 2 \n(l+x)/x is the potential of the host 
halo. Once R v is known, we determine r t as outlined in 
Klypin et al. (1999b). If r t < r max we declare the subhalo 
to be tidally destroyed and remove it from the list of galac- 
tic subhalos. Accreted halos that survive (tdf > h — ta 
and r t > r max ) are assumed to preserve the v c they had 
when they were accreted. 

The resulting velocity function of surviving subhalos 
within the virial radius of the host (~ 200 /i _1 kpc at 
z = 0), averaged over all merger histories, is shown by 
the thin solid line in Figure 1. The error bars reflect the 
measured dispersion among different merger histories. The 
thick straight line is the best fit to the corresponding ve- 
locity function measured in the cosmological N-body sim- 
ulation of KKVP99. The upper dashed lines show the ve- 
locity function if the effects of dynamical friction and/or 
tidal disruption are ignored. The analytic model repro- 
duces the N-body results remarkably well. We have not 
tuned any parameters to obtain this agreement, although 
we noted above that plausible changes in the assumed ini- 
tial circular radius R c could change the analytic predic- 
tion by ~ 50%. The good agreement suggests that our 
analytic model captures the essential physics underlying 
the N-body results. An interesting feature of the model 
is that the subhalos surviving at z = are only a small 
fraction of the halos actually accreted, most of which are 
destroyed by tidal disruption. We discuss implications of 
this satellite destruction in §4. 

2.2. Modeling Observable Satellites 

The second step in our model is to determine which of 
the surviving halos at z = will host observable satel- 
lite galaxies. The key assumption is that after the reion- 
ization rcdshift z ro , gas accretion is suppressed in halos 
with v c < vt- We adopt a threshold of Vt = 30 km s _1 , 
based on the results of Thoul & Weinberg (1996), who 
showed that galaxy formation is suppressed in the pres- 
ence of a photoionizing background for objects smaller 
than ~ 30 km s _1 . This threshold was shown to be in- 
sensitive to the assumed spectral index and amplitude of 
the ionizing background (a similar result was found by 
Quinn et al. 1996). Shapiro et al. (1997; see also Shapiro 
& Raga 2000) and Barkana & Loeb (2000) have suggested 
that very low mass systems (v c < 10 km s _1 ) could lose the 
gas they have already accreted after reionization occurs, 
but we do not consider this possibility here. 

The calculation of the halo velocity function in §2.1 is 
approximate but straightforward, and we have checked its 
validity by comparing to N-body simulations. Determin- 
ing which of these halos are luminous enough to represent 
known dwarf satellites requires more uncertain assump- 
tions about gas cooling and star formation. We adopt a 
simple model that has two free parameters: the reioniza- 
tion redshift z re and the fraction / = M(z rc )/M a of a 



subhalo's mass that must be in place by z ro in order for 
the halo to host an observable galaxy. 

The value of z ro is constrained to z rc ^ 5 by obser- 
vations of high-z quasars (e.g., Songaila et al. 1999) 
and to z ro <J 50 by measurements of small-angle CMB 
anisotropics, assuming typical ranges for the cosmological 
parameters (e.g., Griffiths, Barbosa, & Liddle 1999). The 
value of / is constrained by the requirement that observ- 
able halos have mass-to-light ratios in the range of ob- 
served dwarf satellites. For the subset of (dwarf irregular) 
satellite galaxies with well-determined masses, the mass- 
to-light ratios span the range M/Ly ~ 5 — 30 (Mateo 
1998); dwarf sphcroidals have similar M/Ly, but with a 
broader range and larger observational uncertainties. We 
can estimate M/Ly for model galaxies by assuming that 
they accrete a baryon mass fM a (Qb/Q m ) before z re and 
convert that accreted gas with efficiency e* into a stellar 
population with mass-to-light ratio M*/Ly, obtaining 



M 



= r 1 



Ly 



-1/? 



(3) 



The factor F Q is the fraction of the halo's virial mass M a 
that lies within its final optical radius (which may itself be 
affected by tidal truncation). For the M/Ly values quoted 
above, the optical radius is typically ~ 2kpc (Mateo 1998), 
and representative mass profiles of surviving halos imply 
F' (2kpc) ~ 0.5; however, this factor must be considered 
uncertain at the factor of two level. The value of e» is 
also uncertain because of the uncertain influence of super- 
nova feedback, but by definition e* < 1. Adopting a value 
M*/Ly ~ 0.7 typical for galactic disk stars (Binney & 
Merrifield 1998), n m /D, b ~ 7 (based on Q m = 0.3, h = 0.7, 
and n b h 2 ~ 0.02 from Buries & Tytler [1998]), e* = 0.5, 
and F = 0.5, wc obtain {M/Ly) ~ 5/ _1 . Matching the 
mass-to-light ratios of typical dwarf satellite galaxies then 
implies / ~ 0.3. With the uncertainties described above, 
a range / ~ 0.1 — 0.8 is plausible, and the range in ob- 
served {M/Ly) could reflect in large part the variations in 
/ from galaxy to galaxy. Values of / <J 0.1 would imply 
excessive mass-to-light ratios, unless the factor F Q can be 
much smaller than we have assumed. 

In sum, the two parameters that determine the fraction 
of surviving halos that are observable are z IC and /, with 
plausible values in the range z rc ^5 — 50 and / ~ 0.1 — 0.8. 
For a given subhalo of mass M a and accretion redshift z a , 
we use equation (2.26) of LC93 to probabilistically de- 
termine the redshift Zf when the main progenitor of the 
subhalo was first more massive than Mf — fM a . We asso- 
ciate the subhalo with an observable galactic satellite only 
if z f > z Ic . 

3. RESULTS 

Figure 2 shows results of our model for the specific 
choices of z ro = 8 and / — 0.3. The thin solid line is 
the velocity function of all surviving subhalos at z = 0, 
reproduced from Figure 1. The thick line and shaded re- 
gion shows the average and scatter in the expected num- 
ber of observable satellites (zj > z IC ). The solid triangles 
show the observed satellite galaxies of the Milky Way and 
M31 within radii of 200 /i -1 kpc from the centers of each 
galaxy (note that all results are scaled to a fiducial volume 
l/i~ 3 Mpc 3 ). We see that the theoretical velocity function 
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for visible satellites is consistent with that observed, and 
that the total number of observable systems is ~ 10% of 
the total dark halo abundance. The reason for this differ- 
ence is that most of the halos form after reionization. The 
observable satellites are within halos that formed early, 
corresponding to rare, high peaks in the initial density 
field. 

Other combinations of z rc and / can provide similar re- 
sults because the fraction of mass in place by the epoch 
of reionization is larger if reionization occurs later. We 
find that the following pairs of choices also reproduce the 
observed galactic satellite velocity function: (/, z ro ) = 
(0.1, 12);(0.2, 10);(0.4, 7);(0.5,5). It is encouraging that 
successful parameter choices fall in the range of naturally 
expected values. However, our model fails if reionization 
occurs too early (z rc zl2), since the value of / required 
to reproduce the observed velocity function would imply 
excessively large M/Ly ratios. We conclude that the ob- 
served abundance of satellite galaxies may be explained in 
the CDM scenario if galaxy formation in low-mass halos is 
suppressed after the epoch of reionization, provided that 

4. DISCUSSION 

The suppression of gas accretion by the photoionizing 
background offers an attractive solution to the dwarf satel- 
lite problem. The physical mechanism seems natural, al- 
most inevitable, and requires no fine tuning of the primor- 
dial fluctuation spectrum or properties of the dark mat- 
ter. Although, in principle, a similar explanation could be 
obtained with supernova feedback, reionization naturally 
explains why most subhalos are dark, while the fraction 
~ 10% that accreted a substantial fraction of their mass 
before z IO remain visible today. With feedback, it is not ob- 
vious why any dwarf satellites would survive, and certainly 
not the specific number observed. Photoionization also 
naturally explains why the discrepancy with CDM predic- 
tions appears at v c ~ 30 km s _1 , rather than at a higher or 
lower v c . In contrast, studies of the feedback mechanism 
suggest that it is difficult to achieve the required suppres- 
sion of star formation in halos with v c zl5 km s _1 (Mac 
Low & Ferrara 1998). 

In order to obtain a reasonable M/L ratio for satellite 
galaxies in this scenario, the reionization rcdshift must be 
relatively low, z ro sl2. For higher z ro , it is hard to un- 
derstand how dwarf galaxies with v c < 30 km s" 1 could 
have formed at all, though a blue power spectrum or non- 
Gaussian primordial fluctuations could help, or perhaps 
dwarfs could form by fragmentation within larger proto- 
galaxics. Currently, the reionization rcdshift is constrained 
only within the rather broad range z ro ~ 5 — 50, but it 
might be determined in the future by CMB experiments 
or by spectroscopic studies of luminous high-z objects. 

A clear prediction that distinguishes this model 
from models with suppressed small-scale power or self- 
interacting dark matter is that there should be a large 
number of low-mass subhalos associated with the Milky 
Way and similar galaxies. If we assume that the model 
presented in §2.2 applies in all cases, then the observed 
dwarf satellites should be just the low M/L tail of the 
underlying population. For example, in our fiducial case 
(Figure 2, z rc = 8, / = 0.3), reducing / by a factor of 3 
— corresponding to an increase in the average M/L by a 
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Fig. 2. — Cumulative velocity function of all dark matter 
subhalos surviving at z — (thin solid line) and "observable" 
halos [zi > 2 rc ) (thick solid line with shading), for the specific 
choice of z la — 8 and / = 0.3. The velocity function repre- 
sents the average over 300 merger histories for halos of mass 
M vil (z = 0) = 1.1 x lO 12 ^ 1 Mq. The errorbars and shading 
show the dispersion measured from different merger histories. 
The observed velocity function of satellite galaxies around the 
Milky Way and M31 is shown by triangles. 



factor of 3 — raises the predicted number of galaxies by 
a factor of 6. Reducing / (increasing M/L) by a factor of 
7 raises the predicted number of satellites by a factor of 
10, and accounts for ~ 98% of the dark subhalos. Large 
area, deep imaging surveys may soon be able to reveal faint 
dwarf satellites that lie below current detection limits. 

Dark halo satellites may also be detectable by their grav- 
itational influence. For example, it may be possible to de- 
tect halo substructure via gravitational lensing. Mao & 
Schneider (1998) point out that at least in the case of the 
quadruply imaged QSO B 1422 + 231, substructure may 
be needed in order to account for the observed flux ratios. 
It is also plausible that a large number of dark satellites 
could have destructive effects on disk galaxies (Toth & 
Ostrikcr 1992; Ibata & Razoumov 1998; Weinberg 1998). 
Moore et al. (1999) used a simple calculation based on the 
impulse approximation to estimate the amount of heating 
that the subclumps in their simulations would produce on 
a stellar disk embedded at the halo center. They concluded 
that this type of heating is not problematic for galaxies like 
the Milky Way, but that the presence of galaxies with no 
thick disk, such as NGC 4244, may be problematic. Our 
results indicate that the variation in halo formation histo- 
ries is substantial, suggesting that thin disks could occur 
in recently formed systems with little time for significant 
heating to have occurred, but a detailed investigation of 
this subject is certainly warranted. 

It is tempting to identify some of these subhalos with 
the High Velocity Clouds (HVCs) (Braun & Burton 1999; 
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Blitz ct al. 1998). In our scenario, it would be possible 
to account for some of these objects provided we mod- 
ify the model of §2.2. We currently assume that there is 
no accretion after reionization, and that most of the gas 
that accretes before z rc is converted into stars. However, 
it may be that accretion starts again at low redshift, af- 
ter the level of the UV background drops (e.g., Babul & 
Rees 1992; Kepner, Babul & Spergel 1997), and these late- 
accreting systems might form stars inefficiently and retain 
their gas as HI. If this is the case, HVCs may be associ- 
ated with subhalos that are accreted at late times. We find 
that, on average, ~ 60% of surviving subhalos fall into the 
host halo after z — 1, and ~ 40% fall in after z = 0.5. 

Another prediction is that there should be a diffuse stel- 
lar distribution in the Milky Way halo associated with 
the disruption of many galactic satellites (see Figure 1). 
If we assume that the destroyed subhalos had the same 
stellar content as the surviving halos, then we can esti- 
mate the radial density profile of this component by plac- 
ing the stars from each disrupted halo at the apocenter 
of its orbit, where they would spend most of their time. 
This calculation yields a density profile p*(r) cx r~ a , with 
a = 2.5±0.3, extending from r ~ 10—100 /i~ 1 kpc, which is 
roughly consistent with the distribution of known stellar 
halo populations such as RR Lyrae variables (e.g., Wet- 
terer & McGraw 1996). The normalization of the profile 
is more uncertain, but for the parameters z ro = 8 and 
/ = 0.3, and with the assumption that each halo with 
Zf > z re has a mass in stars of M* = /(f4/f2o)e*M a , we 
find that the stellar mass of the disrupted component is 
- 5 x 10 8 h~ 1 M G . This diffuse distribution could 
make up a large fraction of the stellar halo, perhaps all 
of it. Observationally, it may be difficult to distinguish a 



disrupted population from a stellar halo formed by other 
means, but perhaps phase space substructure may pro- 
vide a useful diagnostic (e.g., Johnston 1998; Helmi et al. 
1999). 

This disrupted population would not be present in mod- 
els with suppressed small scale power or warm dark mat- 
ter. However, it would be expected in the self-interacting 
dark matter scenario. In this case, the distribution would 
probably extend to a larger radius because dark matter 
interactions would disrupt the satellite halos further out. 

There are other problems facing the CDM hypothesis, 
such as the possible disagreement between the predicted 
inner slopes of halo profiles and the rotation curves of 
dwarf and LSB galaxies (Moore et al. 1999; Kravtsov et 
al. 1998; Flores & Primack 1994; Moore 1994). The mech- 
anism proposed here does not solve this problem, though 
more complicated effects of gas dynamics and star forma- 
tion might do so. We have shown that one of the problems 
facing CDM can be resolved by a simple gas dynamical 
mechanism. If this solution is the right one, then the dark 
matter structure of the Milky Way halo resembles a scaled- 
down version of a typical galaxy cluster, but most of the 
low-mass Milky Way subhalos formed too late to accrete 
gas and become observable dwarf galaxies. 
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the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astron- 
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